GLOBAL RESIDUES FOR SPARSE POLYNOMIAL SYSTEMS 



IVAN SOPRUNOV 

Abstract. We consider families of sparse Laurent polynomials fi, ■ ■ ■ ,f n with a 
finite set of common zeroes Zf in the torus T ?l = (C — {0})" . The global residue 
assigns to every Laurent polynomial g the sum of its Grothendieck residues over 
Zf . We present a new symbolic algorithm for computing the global residue as 
a rational function of the coefficients of the fi when the Newton polytopes of 
the fi are full-dimensional. Our results have consequences in sparse polynomial 
interpolation and lattice point enumeration in Minkowski sums of polytopes. 



1. Introduction 

Let /i = ••• = /„ = be a system of Laurent polynomial equations in n variables 
whose Newton polytopes are Ai, . . . , A n . Suppose the solution set Zf in the algebraic 
torus T = (C — {0}) n is finite. This will be true if the coefficients of the fi are generic. 

The global residue assigns to every polynomial g the sum over Zf of Grothendieck 
residues of the meromorphic form 

uj - 9 — A A — 
9 fl ■ ■ ■ fn ti t n 

It is a symmetric function of the solutions and hence depends rationally on the coeffi- 
cients of the system. Global residues were studied by Tsikh, Gelfond and Khovanskii, 
and later by Cattani, Dickenstein and Sturmfels |Ts| IGKh] ICaD| ICaDSj . There are 
numerous applications of global residues ranging from elimination algorithms in com- 
mutative algebra |CaDS| to integral representation formulae in complex analysis |Tsj . 

There have been several approaches to the problem of computing the global residue 
explicitly. Gelfond and Khovanskii considered systems with generically positioned 
Newton polytopes. The latter means that for any linear function £, among the corre- 
sponding extremal faces A^, . . . , A| at least one is a vertex. In this case Gelfond and 
Khovanskii found an explicit formula for the global residue as a Laurent polynomial 
in the coefficients of the system |GKhj . In general, the global residue is not a Laurent 
polynomial and one should not expect a closed formula for it. 

Another approach was taken by Cattani and Dickenstein in the generalized unmixed 
case when all A, are dilates of a single ra -dimensional polytope (see |CaDj ). Their 
algorithm requires computing a certain element in the homogeneous coordinate ring 
of a toric variety whose toric residue equals 1. Then the Codimension Theorem for 
toric varieties |CD| allows global residue computation using Grobner bases techniques. 
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This approach was then extended by D'Andrea and Khetan to a slightly more general 
ample case |DK| . However, it turned out to be a hard (and still open) problem to find 
an element of toric residue 1 for an arbitrary collection of polytopes. 1 

In this paper we present a new algorithm for computing the global residue when 
A l5 ..., A n are arbitrary n -dimensional polytopes (see Section [SJ). The proof of its 
correctness is in the spirit of Cattani and Dickenstein's arguments, but avoids the 
toric residue problem. It relies substantially on the Toric Euler-Jacobi theorem due 
to Khovanskii |Khj . Also in our algorithm we replace Grobner bases computations 
with solving a linear system. This gives an expression of the global residue as a 
quotient of two determinants. The same idea was previously used by D'Andrea and 
Khetan in |DKj . 

There is an intimate relation between the global residue and polynomial interpola- 
tion. In the classical case this was understood already by Kronecker. In Section |3] we 
show how the Toric Euler-Jacobi theorem gives sparse polynomial interpolation. An- 
other consequence of this theorem is a lower bound on the number of interior lattice 
points in the Minkowski sum of n full-dimensional lattice polytopes in K" (Corol- 
lary l^.2j) . Finally, our results give rise to interesting optimization problems about 
Minkowski sums which we discuss in Section |HI 

Acknowledgments. This work was inspired by discussions with Askold Khovanskii 
when I visited him in Toronto. I am grateful to Askold for the valuable input and 
his hospitality. I thank Amit Khetan for numerous discussions on this subject and 
Eduardo Cattani and David Cox for their helpful comments. 

2. Global Residues in C™ and Kronecker interpolation 

We begin with classical results on polynomial interpolation that show how global 
residues come into play. 

Let /i, . . . , f n G C[ii, . . . , t n ] be polynomials of degrees d±, ■ ■ ■ , d n , respectively. 
Let p = ^2di — n denote the critical degree for the fi . We assume that the solution 
set Zf C C n of /i = • • • = f n = consists of a finite number of simple roots. The 
following is a classical problem on polynomial interpolation. 

Problem 1. Given (f> : Zf — > C find a polynomial g G C[ii, . . . , t n ] of degree at most 
p such that g{a) = <p{a) for all a G Zf . 

Kronecker |Krj suggested the following solution to this problem. Since each /j 
vanishes at every a G Zf we can write (non- uniquely) : 

n 

(2- 1 ) fi = y^9ij(tj - aj), a = (ai,...,a n ), g^ G C[ti, . . . , t n }. 

The determinant g a of the polynomial matrix [gij] is a polynomial of degree at most p. 
Notice that the value of g a at t = a equals the value of the Jacobian Jf = det {j^^j 

^For a combinatorial construction of such an element that provides a solution to the problem for 
a wide class of polytopes see |KS| . 
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at t = a, since = 9ij{a) from (|2.1j) . Also g a (of) = for any a' S Zf, a' ^ a. 

Indeed, substituting t = a' into (|2.1j) we get 

n 

= ^5 ij (a')(cij - aj), 
i=i 

which means that the non-zero vector a' — a is in the kernel of the matrix [gij(a')} , 
hence det[gij(a')] = 0. It remains to put 

(a) 



(2-2) g = 



, Jf(a 



9a 



Remark 2.1. When choosing g^ in (|2.1|) one can assume that 

5ij = — h lower degree terms, 

at,- 

where i*^ is the main homogeneous part of /j . Then (|2.2|) becomes 
(2.3) 9 = ( I I <7f + lower degree terms, 

where Jf is the Jacobian of the F{ . 

Definition 2.2. Given g € C[ti, . . . ,t n ] the sum of the local Grothendieck residues 

a G Z / \h---In J 

is called the global residue of g for the system /i = • • • = f n = 0. In the case of 
simple roots of the system we get 

9(a) 



Kf(9) = £ 



aez. 



J f (a) 



Theorem 2.3. (The Euler-Jacobi theorem) Let fi = ■ ■ ■ = /„ = be a generic 
polynomial system with deg(/j) = di . Then for any h of degree less than p = ^ di — n 
the global residue lZf(h) is zero. 

Proof. Consider the function h : Zf — > C, a i— > h(a). According to the previous 
discussion there is a polynomial of the form (|2.3|) which takes the same values on Zf 
as h. In other words, 

h = ( -rr 1 Jf + lower degree terms (mod I), 

y a U jf{a)J 

where we used that the ideal I = (f\, . . . ,f n ) is radical and the roots are simple. 
Comparing the homogeneous parts of degree p in this equation we see that either 
Jp G J, which is equivalent to the Fi having a non-trivial common zero (does not 
happen generically) , or the coefficient of Jf (the global residue of h) is zero. □ 
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3. Global Residues in T n and sparse polynomial interpolation 

Now we will consider sparse polynomial systems and define the global residue in 
this situation. The word "sparse" indicates that instead of fixing the degrees of the 
polynomials we fix their Newton polytopes. The Newton polytope A(/) of a (Laurent) 
polynomial / is defined as the convex hull in Z™ of the exponent vectors of all the 
monomials appearing in / . Note that the Newton polytope of a generic polynomial 
of degree d is a <i-dilate of the standard n-simplex. Such polynomials are usually 
called dense. In what follows when we say generic sparse polynomial we will mean 
that its Newton polytope is fixed and the coefficients are generic. 

Let fx, . . . , f n £ C[tf l , . . . , t^ 1 ] be Laurent polynomials whose Newton polytopes 
Ai, . . . , A n are n -dimensional. We will assume that the solution set Zf C T n of the 
system fx = ■ ■ ■ = f n = is finite. Here T" denotes the n-dimensional algebraic 
torus (C- {0}) n . 

Similar to the affine case define the global residue of a Laurent polynomial g as the 
sum of the local Grothendieck residues 

aeZf 

When the roots of the system are simple we obtain 

9(a) 



. jj{ a y 

where jj = det (tj^zj is the toric Jacobian of the polynomials fx, ■ ■ ■ , f n - 

Notice that Tlj is a linear function on the space of Laurent polynomials and de- 
pends rationally on the coefficients of the fi , since it is symmetric in the roots Zf . 
In Section [21 we will give an algorithm for computing TZj (g) as a rational function of 
the coefficients of the system. 

The following theorem due to A. Khovanskii is a far reaching generalization of the 
Euler-Jacobi theorem. 

Theorem 3.1. |Khj Let /i = ••• = /« = be a generic system of Laurent polyno- 
mials with n-dimensional Newton polytopes Ai, . . . , A n . Let A = Ai + • • • + A n be 
the Minkowski sum. Then 

(1) (Toric Euler-Jacobi) for any Laurent polynomial h whose Newton polytope 
lies in the interior of A the global residue 1Zj(h) is zero. 

(2) (Inversion of Toric Euler-Jacobi) for any <j) : Zf — > C with ^2 a( z Zf 4>(a) = 
there exists a polynomial h whose Newton polytope lies in the interior of A 
such that 4>(a) = h(a) / Jj (a) . 

Let us denote by Sa° the vector space of all Laurent polynomials whose Newton 
polytope lies in the interior A° of A . We have a linear map 

(3.1) A:S A o^d Z f\ h^(^-,aeZ f 

\Jf{a) 
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Then the above theorem says that the image of A is the hyperplane x% = 0} 
in <& z f\ . By Bernstein's theorem the number of solutions \Zf\ is equal to the nor- 
malized mixed volume n! V(Ai, . . . , A n ) of the polytopes [Hj. We thus obtain a lower 
bound on the number of interior lattice points of Minkowski sums. 

Corollary 3.2. Let Ai,...,A n be n- dimensional lattice polytopes in M n and A 
their Minkowski sum. Then the number of lattice points in the interior of A is at 
least n\V(A 1 ,...,A n ) - 1. 

It would be interesting to give a direct geometric proof of this inequality and 
determine all collections Ai, . . . , A n for which the bound is sharp. In the unmixed 
case Ai = • • • = A n = A the inequality becomes 

(3.2) (nA)° n Z n > n! Vol n (A) - 1 

and can be deduced from the Stanley's Positivity theorem for the Ehrhart polynomial 
|5t] . Recently Batyrev and Nill described all possible A which give equality in (|3.2j) 
(see [BN]). Here is a mixed case example which shows that the bound in Corollary 13. 21 
is sharp. 

Example 1. Let T(m) denote the simplex defined as the convex hull in W 1 of n + 1 
points {0, ei, . . . , e n _i, me n } , where is the i-th standard basis vector. 

Consider a collection of n such simplices T(mi), . . . ,T(m n ) with mi < ■ ■ ■ < m n . 
It is not hard to see that their mixed volume equals m n . For example, one can consider 
a generic system with these Newton polytopes and eliminate all but the last variable 
to obtain a univariate polynomial of degree m n . The number of solutions of such a 
system is m n , which is the mixed volume by Bernstein's theorem. 

Also one can see that the number of interior lattice points in T(mi) + • • • + T(m n ) 
is exactly m n — 1. (In fact, these lattice points are precisely the points (1, . . . , 1, k) 
for 1 < k < m n .) 

Corollary 3.3. (Sparse Polynomial Interpolation) Let f± = ■ ■ ■ = f n = be a 

generic system with n-dimensional Newton polytopes A\, . . . , A n and let A be their 
Minkowski sum. Let Zf C T n denote the solution set of the system. Then for any 
function (f> : Zf — > C there is a polynomial g with A(g) C A such that g(a) = 4>(a) . 
Moreover, g can be chosen to be of the form g = h + cjj for some h with A(h) C A° 
and a constant c. 

Proof. Consider a new function tp : Zf — > C given by 

Then the sum of the values of ip over the points of Zf equals zero. Therefore there 
exists h E Sa° such that h(a) = jj(a)tp(a) = <p(a) — jjy Jj(a) for all a G Zf. It 
remains to put g = h + jryJl ■ □ 
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Remark 3.4. Theorem 13.11 is an instance of a more general result. Let /i, •••,/& 
be generic Laurent polynomials with n -dimensional Newton polytopes Ai, . . . , A&, 
for k < n. The set of their common zeroes defines an algebraic variety in T n . 
There is a way to embed T n into a projective toric variety X so that the algebraic 
closure of fi = defines a Cartier divisor D{ on X and the closure Z^ is a complete 
intersection in X . In |Kh2j Khovanskii described the space of top degree holomorphic 
forms on . The special case k = n corresponds to the space of all functions on 
the finite set Z n = Zf whose description we gave in Theorem 13.11 It follows from 
cohomology computation on complete intersections. In particular, for k = n, there is 
an exact sequence of global sections 

n 

(3.3) ► H°(X, 0(D - D { + K)) H°(X, 0(D + K)) -> C 1 ^ 1 C -> 0. 

i=l 

Here D = D\ + • • • + D n , K the canonical divisor, and 0(L) the invertible sheaf 
corresponding to a divisor L . The first non-zero map in (|3.3j) (from the right) is the 
trace map, the second map is the residue map we considered in (|3.1|) . and the third 
one is given by (fx,... ,/„). 

4. Some commutative algebra 

As before consider a system of Laurent polynomial equations fx = ■ ■ ■ = f n = 
whose Newton polytopes Ai, . . . , A n are full-dimensional. For generic coefficients the 
system will have a finite number of simple roots Zf in the torus T n . We concentrate 
on the following problem. 

Problem 2. Given a Laurent polynomial g compute the global residue TZj(g) as a 
rational function of the coefficients of the f . 

We postpone the algorithm to the next section and now formulate our main tool 
for solving the problem. 

Theorem 4.1. Let Aq, . . . , A n be n + 1 full- dimensional lattice polytopes in W 1 and 
assume that Aq contains the origin in its interior. Put 

(4.1) A = A + --- + A n and A (i) = A + • ■ ■ + A;_i + A i+1 + • ■ ■ + A n . 

Then for generic polynomials fi with Newton polytopes Aj , for 1 < i < n, the linear 
map 

n n 

(4.2) © C -> S Ao , (h , ...,h n ,c)^h + Y^ hfi + cJj 

i=0 (l i=X 

is surjective. 

Proof. Let g be in S^ a . By Corollary 13.31 there exists a polynomial ho supported in 
A = A° Q ^ and a constant c such that g(a) = ho(a) + cjj(a) for all a € Zf, i.e. 

the polynomial g — ho — cJj £ S^ a vanishes on Zf . Now the statement follows from 
Theorem E2] b elow . □ 
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The following statement can be considered as the toric version of the classical 
Noether theorem in P n (see |Tsj . section 20.2). 

Theorem 4.2. Let fi,...,f n be generic Laurent polynomials with n- dimensional 
Newton polytopes A\, . . . , A„ . Let h be a Laurent polynomial vanishing on Zf . As- 
sume A(h) lies in the interior of A = Ao + Ai + • • • + A n for some n- dimensional 
polytope Ao . Then h can be written in the form 

h = hifi H h h n f n , with A(hi) C A^, 

where Ar^ as in 

Proof. First we note that the statement remains true when Ao = {0}, i.e. h is 
supported in the interior of A = Ai + ■ ■ ■ + A n . This follows from the exact sequence 
()3.3|) . Indeed, if one considers the toric variety X associated with A then each /j 
defines a semiample divisor Di on X with polytope A,; . It is well-known that for any 
semiample divisor L 

(4.3) H n - dimAL (X,0(L + K)) = S A o l , 

where A^ is the polytope of L (see, for example |Khll IF] ). Thus the first term in 
(13.31) is isomorphic to S<\° © • • • © S&° , where Au-] = Y\^i A,- . Since the sequence 

(1) (n) V / — Jr' J 

is exact and h lies in the kernel of the second map we get the required representation. 

Now assume Ao is n -dimensional. Let X be the toric variety associated with 
A. Let /o be any monomial supported in Ao- Then /o, . . . , f n define n + 1 semi- 
ample divisors Do,... ,D n on X whose polytopes are Aq, • • • , A n . Since f\, . . . , f n 
are generic (and so all their common zeroes lie in T n ) and /o is a monomial, the di- 
visors Do, . . . , D n have empty intersection in X . Then the following twisted Koszul 
complex of sheaves on X is exact (see |CD| IDK] ): 

n n 

0{K) -» Q(Di + K)^ ► 0{D - Di + K) -» 0(D + K) -» 0, 

i=0 i=0 

where D = Do + ■ ■ • + D n and K the canonical divisor on X . The first few term of 
the cohomology sequence are 

n 

► H°{X, 0(D -Di + K)) -> H°(X, 0{D + K)) ^ 0, 

i=0 

where the middle map is given by (/o, . . . , f n ) ■ This, with the help of (|4.3|) . provides 
h = hofo + hifi H \-h n f n , where A{hi) C A^. 

Notice that ho vanishes on Zf and is supported in the interior of Am) = A. Therefore, 
there exist h\ such that 

ho = h'ji + ■■■ + h'j n , with A(h[) c A^ 
by the previous case. It remains to note that A(/i-/o) C A?^ + Aq = A?^ . □ 
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Remark 4.3. Theorem 14 . 1 1 has interpretation in terms of the homogeneous coordi- 
nate ring Sx of the toric variety X associated with A (see [C]). One can homog- 
enize /o, . . . , f n to get elements Fq, . . . , F n £ Sx of big and nef degrees. According 
to the Codimension Theorem of Cox and Dickenstein (see |CD| ) the codimension 
of / = (Fo,...F n ) in critical degree (corresponding to the interior of A) equals 1. 
Then Theorem 14. II savs that the homogenization of the Jacobian J J to the critical 
degree generates the critical degree part of the quotient Sx / 1 ■ 



5. Algorithm for computing the global residue in T n 

Now we will present an algorithm for computing the global residue Tij(g) for 
any Laurent polynomial g assuming that the Newton polytopes of the system are 
full-dimensional. 

Algorithm 1. Let /i = ■■■ = /„ = be a system of Laurent polynomial equations 
with re-dimensional Newton polytopes Ai,...,A n . As before we let A denote the 
Minkowski sum of the polytopes. 

Input: A Laurent polynomial g with Newton polytope A(g). 

Step 1: Choose an re -dimensional polytope Ao with € Aq such that the Minkowski 

sum A = Ao + A contains A(g) in its interior. 
Step 2: Solve the system of linear equations 

n 

(5.1) g = h + J2 h ifi + cJ f 

i=l 

for c, where hi are polynomials with unknown coefficients supported in the 
interior of (see (|4.1(l ). 

Output: The global residue TZj(g) = cnl V(Ai, . . . , A„) . 

Proof. According to Theorem 14. 11 given g with A(g) C A° there exist hi supported 
in A/U and c G C such that g = ho + Y^i=i hifi + C ^J ■ Taking the global residue we 
have 

n 

Kj(g) =Kj(h )+K T f (Y, h ifi) +cnj(jJ) = cn\V(A 1 ,...,A n ), 

i=l 

where the first two terms vanish by Theorem 13. II (1) and the definition of the global 
residue. □ 

Remark 5.1. Notice that we can ignore those terms of J J whose exponents lie in 
the interior of A since their residue is zero by Theorem 13.11 (1), and work with the 
"restriction" of J J to the boundary of A . 

We illustrate the algorithm with a small 2-dimensional example. 
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Example 2. Consider a system of two equations in two unknowns. 

\ fi = aix + a 2 y + a 3 x 2 y 2 , 
\f2 = hx + b 2 xy 2 + b 3 x 2 y 2 . 
The Newton polytopes Ai , A2 and their Minkowski sum A are depicted in Figure ^ 





Figure 1. 

We compute the global residue of g = x 5 y 4 . Let Ao be the triangle with vertices 
(-1, 0) , (0, -1) and (2, 1) . Then the Minkowski sum A = A + A contains A(g) = 
(5,4) in the interior (see Figure EJ). The vector space has dimension 15 and a 




An 





A(<?) 



Figure 2. 



monomial basis 

S&o = (xy, xy 2 ,xy 3 ,x 2 ,x 2 y, x 2 y 2 , x 2 y 3 ,x 3 y, x 3 y 2 , x 3 y 3 ,x 3 y 4 , x 4 y 2 ,x 4 y 3 , x 4 y 4 , x 5 y 4 ). 

Now A( ) = A1+A2 = A, Am = A0+A2 and A( 2 ) = A0+A1 and the corresponding 
vectors spaces are of dimension 4, 6 and 6, respectively. Here are their monomial bases: 



S A O 



I A 22233 3\ o I A A A A 6 A\ 

(x y,x y ,x y ,x y ), S^ = {x,xy,xy ,x y,x y ,x y ), 



2 „2„, 2 ,2 „3,2\ 



"(0) 



as Figure 01 shows. 
We have 



/ 2 2 2 3 2\ 

A ° - \V,x,xy,x y,x y ,x y ), 

(2) 



Jf = -a 2 b\xy - a 2 b 2 xy + 2a\b 2 x y~-2a 2 b 3 x y + 2(ai6 3 - a 3 b\) x y + 2a 3 b 2 x y , 
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Figure 3. 

where we can ignore the third and the fourth terms (see Remark 15. ljl . Now the map 
Q4.2|) written in the above bases has the following 15 x 17 matrix, which we denote 
by A. 



■ 











«2 

















bi 

















-a 2 bi 

















a 2 





















































a 2 











b 2 

















-a 2 b 2 














a i 




















h 

















1 














a i 




















h 

















1 














«i 


a 2 











bo 























1 

















a 2 





h 





b 2 



































"1 

















h 























«3 











a i 








h 








bi 





2(ai& 3 - a 3 bi) 











1 





«3 











a 2 








b 3 


b 2 





























«3 























6 2 





2a 3 b 2 





























a i 

















h 


























«3 

















h 



































«3 

















63 


b 2 





. 


























«3 

















b 3 






It remains to solve the system Ac = b, where c is the vector of unknowns (the 
coefficients of the hi and c) and b is the monomial g = x 5 y 4 written in the basis for 
S Aa , i.e. b = (0, . . . , 0, 1)* . With the help of Maple we obtain 

_1 ai 2 b 2 

4 a 3 (ai b 3 - a 3 bi) 2 ' 

Since the mixed volume (area) of Ai and A 2 equals 4, we conclude that 



_ ai 2 b 2 
1 03(0103-0301) 



6. Some geometry 

The first step of our algorithm is constructing a lattice re -dimensional polytope Ao 
such that the Minkowski sum A = Aq + A contains both A and A (g) in its interior. 
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There are many ways of doing that. For example, one can take Ao to be a sufficiently 
large dilate of A (translated so it contains the origin in the interior). In general, this 
could result in an unnecessarily large dimension of S^ a , which determines the size 
of the linear system (|5.1|) . Therefore, to minimize the size of the linear system one 
would want to solve the following problem. 

Problem 3. Given two lattice polytopes A and A' in R™, dim A = n, find an 
n-dimensional lattice polytope Ao such that A + Ao contains both A and A' in its 
interior and has the smallest possible number of interior lattice points. 

This appears to be a hard optimization problem. Instead we will consider a less 
challenging one. First, since the global residue is linear we can assume that g is a 
monomial, i.e. A(g) is a point. 

Problem 4. Given a convex polytope A and a point u in W 1 find a segment 7 
starting at the origin such that u is contained in the Minkowski sum A + 7 and the 
volume of A + 7 is minimal. 

If 7 = [0, m] , for m € Z n , is such a segment then we can take Ao to be a "narrow" 
polytope with in the interior and m one of the vertices, as in Figure [21 Then the 
volume (and presumably the number of interior lattice points) of A = A + Ao will 
be relatively small. 

We will now show how Problem|I]can be reduced to a linear programming problem. 
Let A C W 1 be an n -dimensional polytope, 7 = [0, v] a segment, v G W 1 . First, notice 
that the volume of A + 7 equals 



where pr^A is the projection of A onto the hyperplane orthogonal to 7, Vol/, the 
k -dimensional volume, and |7| the length of 7. For each facet T C A let nr denote 
the outer normal vector whose length equals the (n — 1 )-dimensional volume of T . 
Then we can write 



rcA 

But the latter is the support function hz of a convex polytope (zonotope) Z , which 
is the Minkowski sum of segments: 



Indeed, hz is the sum of the support functions of the segments. Also it is clear that 
h[- nr ,nr]( v ) = J^. 1 ( tn r,v) = Kn r ,«)|. 

The following figure shows the polytopes A and Z , and the normal fan Yiz of Z . 

Now we get back to Problem^ After translating everything by — u we may assume 
that u is at the origin. Then Problem 0] is equivalent to finding x £ A such that the 
volume of A + [0, —x] is minimal, which by the previous discussion means minimizing 
the support function hz(—x) = hz(x) on A. 



Vol n (A + 7) = Vol n (A) + |7| • Vol n _i(pr 7 A) 




h z(v) = ^2\(n r ,v)\, Z = ^2[-n r ,n r ]. 
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Figure 4. 



We can interpret this geometrically. The function hz is a nonnegative continuous 
function, linear on every cone of the normal fan T*z . Its graph above the polytope A 
is a "convex down" polyhedral set in R n+1 (see Figure [SJ). The set of points with the 
smallest last coordinate is a face of this polyhedral set. The projection of this face to 
A gives the solution to our minimization problem. 




• 

Figure 5. 

Finally, note that the normal fan T,z has a simple description. It is obtained by 
translating all the facet hyperplanes Hy , for T C A , to the origin. 
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